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The nanomechanical single-electron shuttle is a resonant system in which a suspended metallic 
island oscillates between and impacts at two electrodes. This setup holds promise for one-by- 
one electron transport and the establishment of an absolute current standard. While the charge 
transported per oscillation by the nanoscale island will be quantized in the Coulomb blockade regime, 
the frequency of such a shuttle depends sensitively on many parameters, leading to drift and noise. 
Instead of considering the nonlinearities introduced by the impact events as a nuisance, here we 
propose to exploit the resulting nonlinear dynamics to realize a highly precise oscillation frequency 
via synchronization of the shuttle self-oscillations to an external signal. 



Introduction. - Micro- and nanoscale resonators [H [2] 
currently receive much attention both in applied and ba- 
sic research. Their extreme sensitivity and design flexibil- 
ity have enabled, for instance, mass sensing on the single 
ad-atom level [3]-[5], force sensitivity below the attonew- 
ton range [6l [7] , and access to the ultimate quantum me- 
chanical limits of motion [SHTU]. Within this world of 
extremes, nanomechanical electron shuttles have the 
potential to administer current with single-electron accu- 
racy. This ultimate limit of current control would allow 
redefinition of the current standard [l2l - fT4] . closing the 
metrological triangle by joining voltage and resistance 
as exact quantities defined in terms of fundamental con- 
stants. 

A nanomechanical electron shuttle [T5l - f22| consists of 
a nanomechanical resonator carrying a metallic island 
which oscillates between two metallic electrodes (Fig. [T]). 
The island defines a quantum dot which, if operated in 
the Coulomb blockade regime [231 El], can be consistently 
charged with a known number of charge carriers down 
to the single electron level. The advantage of nanome- 
chanical shuttles over other single electron tunneling ap- 
proaches is their inherent suppression of unwanted co- 
tunneling due to tunnel contact with only one electrode 
at a time. Given sufficiently large DC voltage, the shut- 
tle enters into Coulomb attraction-driven self-oscillation, 
where Coulombic forces suffice to drive the charged island 
to impact with and re-charge on subsequent electrodes 
[TBI dni [22], realizing a nanoscale version of Benjamin 
Franklin's "lightning bell" [25] . 

Going beyond the dynamical effects of the alternat- 
ing recharging of the shuttle at the electrodes [211 [23 j 
we model the highly non-linear (and less explored) mo- 
tion in the impacting regime. We show that the impact 
dynamics can be exploited to provide a stable and repro- 




Figure f: (a) A nanomechanical electron shuttle oscillating 
between two electrodes, (b) A typical trajectory of a shut- 
tle driven by an oscillating force and impacting on the two 
electrodes. 

ducible frequency by synchronizing the self-oscillations to 
an external drive. In combination with charge quantiza- 
tion [30] [21] : this will provide the desired reproducibility 
and precision of the shuttle current. 

Model. - The shuttle can be modeled as a nanomechan- 
ical harmonic oscillator with eigenfrequency ujq , effective 
mass rn, and quality factor Q = Wq/F. We note that 
any extra intrinsic nonlinearities, while possibly present, 
do not alter the essential results below, since the main 
nonlinearity is provided by the impacts at the electrodes. 
The synchronization we analyze is a robust phenomenon, 
insensitive in its main aspects to many details of the 
model. The shuttle itself is centered halfway between two 
electrodes that are located at positions ±L and restrain 
its oscillatory motion. Thus, inelastic impact events oc- 
cur at the electrodes, where the velocity reverses sign and 
is reduced by an impact damping parameter ^ < 1. At 
the same time, the shuttle is charged to a value q — ±CV 
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(depending on which electrode it contacts), and after- 
wards is accelerated towards the other electrode by a 
force F — qV/2L (x V"^. A crucial component of our 
discussion is additional external driving, e.g. by a peri- 
odic force. 

The motion between the electrodes is represented by 
the following equation of motion (written in dimension- 
less units): 

Sir) + Q'^z{t) + z(t) = Tcosinr) ± . (1) 

Position z and time r are related to real physical units 
hy X = Lz and r — uj^t. Impacts occur at z = ±1, where 
some kinetic energy will be lost, effectively re-scaling the 
velocity v = z: v' = —yfrjv. The force term in Eq. ([l]) 
carries a sign that depends on which electrode has been 
contacted last ("-I-" for z = —1). 

Our model is described by five dimensionless parame- 
ters: In Eq. ([T]), we introduce (i) the mechanical quality 
factor Q, (ii) the dimensionless amplitude T of the ex- 
ternal driving force, (iii) its external driving frequency 
il, written in units of the oscillator resonance frequency 
Wq, (iv) the impact damping parameter r/ and (v) the 
dimensionless voltage V = V/V^, with CV^ /2 = mwpL^. 
The quantity = {V/V^,)'^ is half the ratio between the 
charging energy and the oscillator's potential energy at 
the impact point. Note that in writing down Eq. ([T]) we 
assume the charge q to be linear in the voltage V . At 
low temperatures, in the Coulomb blockade regime, the 
charge shows discrete plateaus as a function of V. Then 
in Eq. ([l]) and in all of the subsequent formulas and 
figures must simply be replaced by QV / {2muji^L?) . Fig- 
ure [l] displays a typical trajectory with impact at the 
electrodes, under simultaneous driving. 

All the physical quantities of interest may be expressed 
by the dimensionless parameters and a few dimensional 
quantities. For example, the average shuttle current can 
be obtained from the transported charge CV and the 
frequency: / = CVuJof~^. Here is defined as the 
averaged inverse (dimensionless) time for a one-way trip 
between the electrodes. 

Self-oscillations without external driving. - When a 
sufficiently large dc voltage is applied to the electrodes, 
the shuttle can execute self-oscillations even in the ab- 
sence of external resonant driving This happens 
when the acceleration by the electric field overcomes 
the frictional losses which are mostly due to the impact 
damping. While the existence of this regime is well- 
known, we will discuss it here, since its properties are 
crucial for our subsequent analysis. In this regime, the 
shuttle shows regular motion, where the velocity before 
impact is always the same. We first assume intrinsic os- 
cillator losses to be absent (i.e. Q — oo, but r] ^ 0). 
The energy lost upon impact is obtained via the relation 
v' = yjr]v between the speeds before (v) and after (w') 
impact: (v^ — w'^)/2 = v^iX — ri)/2. This must equal the 
energy gained by acceleration in the electric field, 2V^ in 



our units. Equating the two yields the velocity before 
impact: 

^ 4VV(1 - . (2) 

Now we can obtain the one-way travel time tq. We as- 
sume the last impact was at z(0) = —1, with i(0+) = u'. 
Based on the solution z(r) = + w'sin(r) — (1 + 
V^)cos(t), we demand z{to) = +1. This yields tq via 
sin To = - 4AC]/2A , with A = v'^ + {l+V^)^, 

B = 2v'{V^ - 1), and C = -4V^. 

We note the following important physical features. De- 
spite the impact damping, self-oscillations are possible 
down to the lowest voltages for Q = oo, because the 
grazing impact of the shuttle minimizes the loss: im- 
pact velocity v for V — )• 0. In this limit, the one- 
way travel time is just half the intrinsic oscillator period 
(tq = tt), i.e. the impact happens exactly at the turning 
point of the oscillatory motion, with the shuttle barely 
touching the electrode. Consequently, tq then turns out 
to be independent of the impact damping parameter rj. 
For small voltages, we obtain tt — r « 2V(l-|-y^)/-y/l — rj. 
This relation represents the decrease of the travel time 
from its zero-voltage value, to first order in V, and it 
could serve to extract the impact damping rj. At high 
voltages (V 3> 1), one would observe / cx V^, as higher 
accelerations lead to shorter travel times. However, cur- 
rent experimental parameters 22 indicate that typically 
^ 10~2, which suggests that the high- voltage regime 
is not reached. 

We now reconsider the intrinsic oscillator damping 
(thus Q < oo), where a minimum voltage has to be ap- 
plied for self-oscillations [TTl [53] in order to overcome 
the frictional losses. This voltage can be obtained by 
demanding the total loss (both by impact damping and 
mechanical friction) during one half-cycle to equal the en- 
ergy gain by electrostatic acceleration. This then defines 
the threshold dc voltage for the shuttle's self-oscillations 
in the absence of driving: 

vL-^Q-^ (3) 

Numerical calculations for the experimentally relevant 
parameters suggest that apart from the appearance of 
this threshold voltage there is no appreciable modifica- 
tion of the shuttle dynamics above threshold as long as 
Q is large ^2] . 

General nonlinear map in the presence of driving. - 
Much more complex nonlinear behavior can be observed 
in the presence of an additional oscillatory drive, where 
several different attractors or even chaotic motion at 
higher drive strengths may arise. 

Suppose an impact occurs at z = —1, with an initial 
velocity u„. Once we also specify the phase tpn = ilr'") 
of the external oscillating force at the time t'"^ of impact. 
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Figure 2: Nonlinear map for the impact dynamics of an elec- 
tron shuttle. Phase of the external drive and impact veloc- 
ity V axe displayed at each impact event, (a) The evolution 
maps (}Pn,Vn) into i;„+i) upon the next impact at the 

same electrode. This map completely encodes the dynamics, 
(b) Fixed-points, both attractive and repulsive, arise in the 
[ip, «)-plane for suitable parameters. These are intersections 
of the curves where Vn+i = or ipn+\ = fn (shown here). 
Small parameter changes usually do not change the number 
and type of fixed points, which makes synchronization robust. 



these two values (?;„, (/?„) completely determine the sub- 
sequent shuttle evolution. Evolving the shuttle towards 
the next impact at z = — 1 (usually with an intermediate 
impact at z = 1), we find new values (pn+i and Vn+i, 
which are determined by a unique mapping (cf. Fig. [2]): 



Vn+1 = v'{(pn,Vn) 



(4) 

(5) 



It is known that such two-dimensional nonlinear map- 
pings can generate both complex attractors with a pe- 
riod larger than one cycle and chaotic dynamics. Below, 
we will be interested in synchronization, which (in the 
simplest case) is described by a period-one fixed point 
of the type (p = ip'{ip,v), v = v'{(p,v). Note that in 
the absence of an external drive, there would be only a 
one-dimensional relation w„+i = v'{vn) which goes into a 
simple fixed-point. That is the regular shuttling already 
discussed above. 

Synchronization to an external drive. - Although the 
map of Eqs. ([4j [s]) can generate a wealth of phenomena, 
our analysis below focuses on the case of greatest po- 
tential impact for applications. This is synchronization 
to an external drive (injection locking), which may be 
exploited to lock the shuttle dynamics to a very precise 
external frequency source even for a rather weak drive. 
In this way the nonlinear dynamics can be turned from 
a complication into a useful tool. 

At small external driving strength, there is no simple 
fixed-point. This is because the phase just proceeds ac- 
cording to ipn+i ~ + 2flTQ, with tq the one-way travel 
time in the absence of driving. Generally, 2toQ is not an 
integer multiple of 27r, and so there is no definite phase 
relation between the external driving and the shuttling. 
The combination of the shuttle's impact dynamics and 
the weak incommensurate driving makes the shuttle mo- 
tion lose its exact periodicity. As a result, the shuttle's 
round-trip time 2to effectively changes slightly from step 
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Figure 3: Synchronization of an impacting electron shuttle 
to an external drive, (a) The phase increment ipn+i — ifin 
of the external drive between two impacting events is shown 
as a color-code. As time progresses (vertical axis, measured 
in shuttle cycles), the phase increment either displays a slow 
beating pattern, or it may lock to a value of 2it, signalling syn- 
chronization. This happens inside a certain interval, as the 
dc shuttle voltage (or any other parameter) is varied, shifting 
the shuttle's eigenfrequency. [r? = 0.2, Q = 10^, J" = 10"^, 
57 = 1.054 ] (b) The effective, time-averaged shuttle frequency 
Ld — n/f is plotted vs. an external parameter (the voltage). 
Inside the synchronization interval, the shuttle becomes ro- 
bust against noise (brown curve) - the curve shown here was 
obtained for fluctuations Ail = 0.02 in the driving frequency 
(see main text) and is evaluated after 6000 impact events. 
The fiuctuations decrease for longer averaging times. 



to step. Fig. [3^ depicts the phase increment (pn+i — y^n as 
a function of time and shows long-term beating patterns. 
However, when changing the shuttle's parameters (e.g. 
the voltage), the shuttle's intrinsic period shifts until it 
becomes almost commensurate with the external drive 
frequency. Then synchronization may set in (Fig. [S]). 

Mathematically, this corresponds to the generation of 
a stable period-one fixed point. We now examine the ap- 
pearance of this phenomenon both numerically and an- 
alytically. Numerical observations indicate that the sys- 
tem, once started anywhere in the {ip, u)-plane, quickly 
relaxes to a Id manifold that can be described by the 
phase coordinate alone (Fig. |4^). Therefore, we focus on 
the phase map and expand it as 



ipn+i = fn + 2fiTo -I- K sin((p„ 



p* 



(6) 



Here the values of K and (/3* depend on the detailed mi- 
croscopic parameters, with the coupling K growing from 
zero upon increasing the drive, and the omitted terms 
are higher harmonics. Let us denote the deviation of the 
drive frequency from the undriven shuttling frequency as 
50, — VI — tt/tq. If this is small enough, then we obtain a 
fixed point [fn+i — H^n = with some phase lag (p be- 
tween shuttling and drive, where 2x0(50 = Ksva{(p — Lp.^). 
This requires |2riTo — 2tt\ = \25VItq\ < K, which defines 
the synchronization interval, i.e. the permissible devia- 
tion Sfl. Outside that interval, the phase drifts across 
the full range [0, 27r], with a varying velocity. Near syn- 
chronization, we can assume small phase increments, and 
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Figure 4: Phase dynamics in an impacting electron shuttle, 
(a) After only a few impact events, the system dynamics 
quickly settles onto a Id manifold, permitting a theory for 
the phase dynamics alone, (b) The resulting Adler equation 
corresponds to a particle sliding in a washboard potential, 
(c) The predicted synchronization interval (bounded by the 
contour lines shown here, partially hatched) as a function of 
drive frequency, plotted using our approximate analytical pre- 
diction for the effective Adler coupling constant of a driven 
impact shuttle, Eq. [r] = 0.7] 



the phase map of Eq. ^ can be written as a differen- 
tial equation by replacing ipn+i — iy9„ — 27r ss 2Tod(p/dT, 
introducing the time-dependent phase shift ^{t). This 
yields a generic equation proposed by Adler to treat 
phase- locking phenomena [281 El] . Integrating the Adler 
equation exactly one finds periodic dynamics for the 
phase increment in each cycle, with a period diverging as 
|2ro<5ri — K\^^^'^ near synchronization. All of this can be 
observed in the direct simulations (Fig. [s]) of the nonlin- 
ear shuttle dynamics under drive and dc voltage. Outside 
the synchronization interval we see a periodic modulation 
with the period diverging near the onset of synchroniza- 
tion, as predicted. 

While Adler-type equations appear generically in syn- 
chronization physics, they becomes predictive only when 
we are able to connect the coupling K to the microscopic 
parameters of the specific model. To this end we repeat 
the calculation of the one-way travel time for the driven 
oscillator, obtaining the correction St to the travel time 
To in the leading order of the driving J^. This also in- 
volves finding the change in the impact velocity by a 
new self-consistency equation, where we exploit the fact 
mentioned above that the velocity quickly relaxes to the 
limit cycle while the phase dynamics are much slower. We 
omit the rather lengthy calculation and just state that St 
becomes a harmonic function of ip. The corresponding 
correction to the phase increment is then Sip — 2nST{(p). 
From this, we obtain the Adler coupling constant K: 



K 



T_ 

2^ 



MU - vqU 



U 



(7) 



1 - ^(cos(to) - Mg) 
Here v = 2V/ y^l — rj is the impact velocity without drive 



(see above), and we abbreviated U = c'-'^^tq — sin(ro), 
U = 2e*'^"— cos(to), g = sin(ro)/w, and M = — w'sin(To) + 
(1 + V^)cos(ro). Note that Eq. ^ is given under the 
simplifying but realistic assumptions of large Q 3> 1 and 
near-resonant driving « 1. We find that Eq. ^ accu- 
rately matches the numerical results for sufficiently low 
drive 

The synchronization interval is given by \Sfl\ < 
K/{2to), growing linearly with drive strength. As can 
be seen in Fig. [4j;, the synchronization interval grows 
for lower voltages (while shifting towards lower frequen- 
cies). At very small voltages, modifications occur due 
to intrinsic damping (finite Q), which destroys the self- 
oscillations. 

Stability against noise. - The existence of an extended 
synchronization interval implies an increased robustness 
of the shuttling against noise. Within this interval, the 
Adler equation corresponds to a phase particle trapped 
in a local minimum of a washboard potential (Fig. |4j3). 
There the effective shuttle frequency is completely in- 
sensitive to a slow drift of system parameters (voltage, 
effective capacitance, damping, etc.). In contrast, time- 
dependent noise may give rise to phase slips by lifting 
the particle over the barrier. However, this process is ex- 
ponentially suppressed when moving towards the middle 
of the synchronization plateau [30] ■ We have confirmed 
this by direct simulation, for both (i) a noisy driving 
frequency fl (reset at each impact to a new Gaussian 
random value) and (ii) a noisy impact parameter rj. Our 
simulations do not show any noticeable effect of a noisy 
T] on the synchronization range, even for Ar]/{r]) 0(1)- 
A leading-order perturbative calculation of the noise- 
induced correction to the one way travel time confirms 
the ineffectiveness of impact noise. For a noisy we 
observe phase slips which drive the effective shuttle fre- 
quency u! away from (S7). They are numerous for volt- 
ages at the edges of the former synchronization interval. 
However, in the center of that interval (e.g. V=0.053 in 
Fig. |3]), for a noise value of Ail/fl = 0.02, we observe 
no phase slips within 300,000 iterations, resulting in a 
frequency imprecision of at most 10~^ in this example. 
The remaining fluctuations of the effective time-averaged 
frequency uj decay with time as l/^/r, as expected. Devi- 
ations of u} from the drive fl are exponentially suppressed 
in the middle of the interval. 

Realization. - Figure [ija) depicts an experimental 
realization that strongly suppresses cotunneling events. 
Here, doubly-clamped, tensile-stressed silicon nitride 
string resonators with linear quality factors of several 
100,000 [311 132] can be be efficiently and controUably 
driven to sufficiently high (impacting) amplitudes of sev- 
eral lOnm using dielectric gradient field actuation [33] 
(acoustic [17] or capacitive drive [151 E] can be used 
in similar contexts). Thus, during tunneling, the is- 
land is always tens of nanometers away from at least 
one electrode. The requirement that the charging en- 
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ergy dominates fc^T can be satisfied with an island- 
electrode capacitance in the attofarad range, imposing 
experimentally accessible typical island cross sections of 
tens of nanometers. With these prerequisites within 
reach [T71 [33], the ideas presented here would enable 
quantitatively reliable current measurements with a me- 
chanical shuttle. 

Conclusions. - In this Letter we have analyzed 
nanomechanical single electron shuttles in the impact 
regime. We predict that synchronization to an exter- 
nal driving frequency can occur and can be exploited to 
achieve a very precise shuttle current. This has imme- 
diate applications in the active experimental research on 
mechanical electron shuttles. In addition, the framework 
presented here can form the basis for discussing higher- 
order fixed-points and chaotic motion in shuttles and sim- 
ilar impacting systems. 
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